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A generalized optical potential for elastic scattering induced by light nuclei is calculated within 
the Feshbach projection operator method. The model explicitly takes into account the contribution 
of the projectile break-up continuum treated within a microscopic few-cluster model. In this work 
we formulate the model, deriving an explicit expression for the optical potential, and show ability 
of the model applying it to deuteron elastic scattering. 



PACS numbers: 24.10.-i, 24.10.H, 25.45.-z, 25.60.Bx 



I. INTRODUCTION 

During the last years reactions with light weakly bound 
nuclei have been of increased interest from the experi- 
mental and theoretical point of views and the progress 
in the investigation of these nuclei has been impressive 
[J, [2] • This progress is conditioned by the important ef- 
forts devoted to investigate reaction mechanisms and by 
new techniques enabling the production of exotic nuclei. 
In spite of that we are still far from a clear understanding 
of the unusual structure of exotic nuclei and of the reac- 
tion mechanisms induced by these nuclei. This arises 
from the experimental difficulties (also by low intensi- 
ties of available beams) and by the difficulties arising in 
the description of a few-body nuclear dynamics taking 
place under conditions of a strong coupling of all reac- 
tion channels with the break-up channel of weakly bound 
projectiles. 

Generally it is assumed that in light exotic nuclei the 
nucleons tend to group into clusters, whose relative mo- 
tion mainly defines the properties of these nuclei. This 
assumption leads to great advantages for models employ- 
ing the cluster concept both for the structure and reac- 
tions involvin g lig ht exotic nuclei [!, 0, H, H, 0, fl & 03 
HH [l2l [H, H3.1l5|. In the case of nuclear reaction study 
the coupled-channel (CC) formalism is the one of con- 
sistent and efficient approaches [H, E3, EH Ell H3] • The 
continuum-discretized reduction of this method (CDCC) 
allows to study the interplay between elastic and break- 
up channels in reactions involving two-cluster nuclei 
EH El, El El El El EH- However, the application of 
the CC approach to reactions with few-body projectiles 
meets significant computational difficulties if a realistic 
few-body wave function is used. Consequently very few 
works have been done [H [H| . 

The generalized optical model (GOM) of H. Feshbach 
[24| is an alternative approach to the problem. Models 
based on the Feshbach theory are extensively used for 
the study of coupling effects on different reaction chan- 
nels. It is worth to mention, in particular, the studies 



of the role of the deuteron break-up in its elastic scat- 
tering with heavy ions [H [H [27], HH, the influence 
of collective excitations on heavy ion elastic scattering 
[H [H EH, EI [3!, [li|, the interplay between break- 
up and complete fusion channels in weakly bound nucleus 
reactions [3a|. However the application of the Feshbach 
method was hampered in the past by the complexity of 
the formulation and of computational burden. As a re- 
sult, even in studies of the deuteron break-up, the appli- 
cations of this method have been done with many sim- 
plifications. 

In contrast with the CDCC approach, the GOM al- 
lows one to avoid the simplifying discretization of the 
continuum spectrum. The calculation within the GOM 
is faster. On the other hand, the CDCC method pro- 
vides phase-shifts in all reaction channels, rather treats 
the continuum-to-continuum coupling, and does not re- 
quire an approximation for the Green function, while the 
GOM does. Therefore in the case of a two-cluster pro- 
jectile the CDCC approach is somewhat more efficient. 
However, in the case of a few-body projectile the appli- 
cation of the CDCC method becomes difficult (the large 
size of the coupling matrix, the complicate procedure of 
matrix element calculation and e.g.), while the GOM re- 
mains more feasible. 

Our main goal is to apply the GOM to study the elastic 
scattering of light weakly bound nuclei using a realistic 
few-body model description of their internal structure. 
As a test of the approach, in the present work, the GOM 
is applied to the deuteron elastic scattering from heavy 
targets at intermediate energies in order to draw conclu- 
sions about its applicability. Applications of the model 
to reactions with few-cluster weakly bound nuclei (such 
as 6 He) will be done in a subsequent publication. 

The method proposed earlier [H [27| is extended here 
in order to avoid simplifying assumptions. In [2H [27j 
the authors (i) neglected the spins of the particles, (ii) 
considered coupling with s-wave continuum only, (iii) ne- 
glected or treated the Coulomb forces in an approximate 
way, and (iv) used the free-particle Green function in- 
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stead of the total one. Our approach goes beyond these 
assumptions. Within the method a structureless target 
nucleus interacts with a projectile treated as a system of 
few bound clusters. The bound and continuum states of 
the projectile are described in the framework of the mi- 
croscopic cluster model and used to construct the Fesh- 
bach projection operators. We derive an explicit expres- 
sion for the optical potential which takes into account 
explicitly the coupling with projectile break-up channels. 
We show also the importance of an accurate treatment 
of the Green function appearing in the dynamical polar- 
ization potential in the case of light targets. 



II. MODEL 

We consider the scattering of a weakly bound projec- 
tile by a structureless target. The projectile is treated as 
a bound few-cluster system. The corresponding Hamil- 
tonian has the following form 



(1) 



where Tr is the kinetic energy operator of the projectile- 
target relative motion, TL^ = t £ + i>£ is the Hamiltonian 
describing the projectile internal structure, and £ denotes 
an appropriate set of internal coordinates. The interac- 
tion potential Vr^ is a sum of effective (non-hermitian) 
cluster-target potentials, which is obtained from a fit of 
elastic scattering data. By employing these cluster-target 
interactions we implicitly take into account the internal 
properties of the clusters and target. 

The total scattering wave function | < E > )r. i ^ satisfies 
the Schrodinger equation with the Hamiltonian (fTJ) and 
eigenvalue E = E p + eq, where E p — H 2 p 2 /2m is 
the projectile-target relative energy, m is the reduced 
projectile-target mass, while £o is the projectile ground 
state energy. The ground \eo : jo)$, (here jo is a projectile 
total angular momentum) and continuum states of 
the projectile form together the spectrum of the Hamil- 
tonian H^, since the weakly bound projectile is supposed 
to have only one bound state. The Feshbach projection 
operators and [24[ are constructed as follows 
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The operator extracts the elastic component of the 

total wave function -P^j'I^R^ = | v I / p)r,j, which satisfies 
the Schrodinger equation 

(Tr + [PVk*P + PV R< ,Q—^QV^P])\^) ntC 



(P-£ )|*p>R,C 



(2) 



where the expression in the square brackets is the gener- 
alized optical potential. The first term U^ 1 ' = PVr.^P is 
the local cluster-target interactions folded over the pro- 
jectile ground state. The second term (we will refer to it 



as U^) is the non-local dynamical polarization potential 
(DPP), describing the coupling of elastic and non-elastic 
channels. 



A. Optical potential for iV-cluster projectiles 

Following the usual technique [37J the elastic compo- 
nent of the many-body wave function may be expanded 
in partial waves as 
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^1+1 i> JLjo (p,R)Q L Z o (Z,n R ), (3) 



where p is supposed to be parallel to the z-axis. 
ipjLj (p,R) is a partial wave function describing the 
projectile-target relative motion with total (J) and or- 
bital (L) angular momenta, and with asymptotic 



i> v (p,R^oo)^F L {pR) 
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H^ipR), (4) 



where S„ = 



„2i<5„ 



H 



(+) 



e"""" is a scattering S-matrix element, 
— Gl + iF L , while F L and Gl are the regular 
and irregular Coulomb functions. In Eq.Q we use the 
notation 6j^*}(£, Qr) for the spin-angle wave function re- 
sulting from the L and jo vector coupling 
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(5) 



where 0j o7o (£) y]° mo (^i) is a partial component of the 
projectile ground state wave function. 

We substitute expansion in the Schrodinger 

equation @. Multiplying the resulting equation by 

®j'M'( r i ^r) fr° m the right, and integrating over r and 
fl r, one obtains a set of coupled Schrodinger equations 
for the partial wave functions 

+R J R'U^l,(R,R')ip v ,{R')dR') =0, (6) 



where index v denotes the set of quantum numbers 
{JLjo}. The sets v and v' differ in orbital momentum 
L and V only, and are used here just for the sake of 
simplicity. Note that the non-diagonality is very weak in 
the case of reactions considered here, and will be finally 
neglected (see below). Thus the set of equations © be- 
comes uncoupled. In the numerical procedure we reduce 
each integro-differential equation ([6]) to a set of linear al- 
gebraic equations applying the finite difference and the 
Simpson methods to approximate the second derivative 
and the integral, respectively. 
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1. Cluster-folding potential 

The projectile-target interaction Vr^ is chosen as a 
sum of effective complex cluster-target potentials. Each 
potential is supposed to be a sum of Coulomb, nuclear 
and spin-orbit terms. The cluster-folding potential 
can be defined as the integral 

u { lI(R) = I d^dn R Q^ L M jo ^,n R )v(-R,oe L jZ(^R)- 

(7) 

The cluster model has been used intensively in stud- 
ies of reactions involving light nuclei. Thus an explicit 
expression for the cluster-folding potential in some par- 
ticular cases can be found elsewhere [H, HH . Therefore 
we give here its final expression only. Using the Fourier- 
Bessel transform of the potential Vr,£, and performing 
the integration over angles and the summation over the 
momenta projections, we obtain the cluster-folding po- 
tential in the following form 

u { L l(R) = 

^J2(-) X L'\ 2 C^ XQ W(j JXL; L'j )F joXjo (R), (8) 

A 

where a = (2a + l) 1 / 2 , C^ ab g is a Clebsch-Gordan co- 
efficient, W(abcd; ef) is a Racah coefficient. The radial 
factor in Eq.© reads 

N J 

Fj„\ 3o (R) = ^2Y1 / V 2 j*( ( lR)vn{q)(jolo\\C n ,x\\jojo)dq, 

™ =1 7o7o o 

(9) 

where N is a number of clusters, v n (q) is a Fourier trans- 
form of the cluster-target potential, and operator C n> \^ 
is defined as 

C nM = X- 1 V^jx(qr n )Y Xfi (n n ). (10) 

The reduced matrix element arising in Eq. ([9]) is a partial 
component of elastic form-factor of n-th projectile cluster 
(i.e. the Fourier transform of its spatial distribution). 

The contribution of the spin-orbit interactions v[ SOS> (r) 
to the folding potential can be calculated in the same 
way as the central one. However both phenomenologi- 
cal and theoretical analyses [4(J EH show that the spin- 
orbit interaction plays a minor role even for deuteron 
elastic scattering and becomes almost negligible in the 
case of heavy ion elastic scattering 41]. It was shown 
[HI that the nucleon-nucleus spin-orbit interaction in the 
first-order perturbation also gives only a spin-orbit term 
for the potential of s-wave projectile (such as deuteron 
or 6 Li). It allows us to consider the folding spin-orbit po- 



G(z) = [z-Ht-f R - Va.e)" 1 . We write the total 
Green operator in the form of the Born series as follows 



G = g + g{v R ^ -v)g 



(ii) 



tential within the approximation proposed in [251 l38l. |41| ] . 



2. Dynamical polarization potential 

The calculation of the polarization poten- 
tial UW(R,R') requires the definition of ma- 
trix elements of the many-body Green operator 



where the operator Q = (z — Tr — TL^ — V) -1 may 
be factorized if the potential V depends on the relative 
projectile-target coordinates R only. The appropriate 
choice of the potential V(R) is the cluster folding poten- 
tial U^ X \R) obtained above. Note, that in this case Hr 
is a nonhcrmitian operator because the folding potential 
U^(R) = U(R) + iW(R) is complex. 

For the total Green operator we use the approxima- 
tion G » Q. By this we neglect the transfer channels, 
whose contribution is included implicitly through the ef- 
fective cluster-target interactions. Omitting the second 
term we neglect also the contribution from the multi- 
step processes like continuum-to-continuum excitations. 
The role of continuum-to-continuum coupling has been 
investigated, in particular, in the paper of Sakuragi et al. 
[ItJ within the CDCC approach. Significant contribu- 
tion of the continuum-to-continuum coupling was shown. 
However, one should point out differences in the formula- 
tions of the models in our work and in Ref. [13] • The au- 
thors of Ref.[l7| defined the projectile-target interaction 
by using the M3Y-type effective nucleon-nucleon interac- 
tion with complex normalizing coefficient and applied a 
double-folding procedure for the calculation of coupling 
matrix elements. The DPP obtained in Ref. [13] has an 
additional large repulsion of the real part in the periph- 
eral region (Re ~ 20 MeV) and an almost negligible 
additional contribution to the imaginary part. On the 
other hand, in Ref.(42j the DPP for the similar reaction 
was calculated within the adiabatic breakup model [43| . 
which is q uite close to our approach. The DPP found 
in Ref.[42| provides weak real (Re < 3 MeV) and 
strong absorptive (7m ~ —10 MeV) contributions. 
Our results are similar. 

In Ref.fl5| the continuum-to-continuum coupling was 
also neglected. The authors concluded that this approx- 
imation should be valid if the elastic component of the 
total wave function is larger than break-up ones. Obvi- 
ously this condition becomes stronger with higher colli- 
sion energies. 

It is difficult to define an explicit applicability condi- 
tion for our approximation. We need to compare the 
collision time t co u with the excitation-deexcitation time 
which is unknown. Instead, we may use the time tint 
associated with the cluster-cluster relative motion inside 
the projectile. This gives us the validity criteria in the 
following form r = t co u/t; mt < 1. For all reactions con- 
sidered in this work this condition is fulfilled. 

Since the total Hamiltonian Hr^ ~ Tr + Vr + 
is separable, one constructs a basis as the direct prod- 
uct of the Hr and fi^ bases. The Green operator Q 
may be then decomposed over the few-body partial states 
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\EpEkJMLjj) which reads 



1/2 



(£,R\E p e K JMLn) =C s i L e^ (f^) 

x^j Lj (p,R)E L J ^\ K ,C,Q R ), (12) 

where e K is an energy of the cluster-cluster relative mo- 
tions, and Q is a phase-volume quantity, which depends 
on the model used for the description of the projectile 
structure. The spin-angular wave function 



JM 



has the same structure as function ©jm(£,^r) except 
the sum over 7. 

The dynamical polarization potential is written as 

00 

ul 2 M',R) = EI\c s \ 2 

Vi 

xVt„. ( K , R') gif (E - e K - R', R) Y ViU . («, R)de K> (13) 
where the matrix elements read 

^vui{n,R) = S JJi 5 MMi 

(14) 

and the quantum number sets v — {JMLj } and ^j_= 
{ Jj-MjijjiTj}. The partial two-body Green function [34| 
in Eq. (T5H is 



gi+\E p ;R,R') 



2m 1 y>„.(p, fl < )fej+ ) (p, iZ>) 



where notations i?< and i?> refer to the smallest and 
largest of coordinates R and R'. The wave func- 
tion (p, i?) = 93^ (p, R) + iip Ui (p, R) is a combina- 
tion of two linear-independent solutions of the two-body 
Schrodinger equation with potential V(R) = f/W(R). 
The regular wave function tp Vi (p, R) has the boundary 
conditions Q , while the irregular solution <p Vi (p, R) has 
the following asymptotic form 



<p Vi (p,R^0) 

in ) 



(pR)- 1 , 
G L { P R n 



1 



H^\pR r , 



(15) 



The correct calculation of the Green function is impor- 
tant for the GOP calculation, since it defines the radial 
dependance of the polarization potential. Therefore the 
free-particle Green function, which was used in [25|, [2tJ , 
provides a crude approximation and could be used for a 
qualitative analysis only. Note also, that the local plane 
wave approximation for the Green function [3] gives sim- 
ple and fast calculation method, which is however less 
accurate than Eq.(TT5]), but still applicable. 



The matrix elements Y VVi {k, R) have the same struc- 
ture as the integral ((7]) except the difference between 
functions and 5. Following the scheme used in the 
case of the cluster-folding potential, one gets the polar- 
ization potential as 



U$ L {R',R) = ^EE(-) V+A??i '^ C ™'oCAo WijoJX'L^L'ji) W&JXL; L lJo ) 



Lai A' A 



^2 J IQI 2 • ;r O'o7o)^'fe'v«)( £ «)- R ') 9j + l]( e ~ £ K ;R',R) F(jm)Mjoio)( £ Ki R ) de *- 



7 7i7o 



(16) 



where the functions T(e K: R) read 

Y 00 

E JV 3x(qR) v„{q) (jo^o\\C n ,x\\s K ,jai)dq. (17) 

n=l 

The reduced matrix element in Eq. (|17[) is a component 
of the transition form-factor, which depends on the pro- 
jectile structure only and may be calculated once. The 
integrand in (|17[) is an oscillating and decreasing func- 



tion of q, and the integration can be done quite easily 
with a truncation at q m ax ~ 6 fm" 1 . The resulting func- 
tion T{e K , R) describes the transition probability and de- 
creases with increasing e K . Therefore an integration in 
Eq. irTB)) may be performed up to some appropriate pro- 
jectile excitation energy £/., which is about 40 MeV in 
the deuteron case [iH ]. 

The polarization potential (fT5|) is non-diagonal on L 
and L 1 indexes. Consequently, one needs to solve numer- 
ically the system of coupled integro-differential equations 
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([6]). It can be done by an iteration procedure using a 
solution of the uncoupled system as an initial approxi- 
mation. However, we avoid this complicated procedure 
here. The potential (|16p has a more simple form in the 
case of the s-wave projectile (like deuteron or 6 Li). It 
is still non-diagonal [V = L, L ± 2, . . . , L ± 2 jo), and 
this non-diagonality is arising from the non-central part 
of the cluster-cluster interaction. The calculations show 
that the non-diagonal terms are about 10 2 times smaller 
than the diagonal one, and therefore can be neglected. 

3. Two-cluster projectile's form-factor 

Let us define the bound and continuum state wave 
function treating the projectile as a two cluster system 
(e.g. the deuteron d=p+n). Vector r = n — r 2 is an ap- 
propriate choice of the coordinate £ for the description of 
the cluster-cluster dynamics. A wave function describ- 
ing bound and scattering states of a two-body system 
e, jmls) (7 = Is) at a relative energy e can be expressed 
as 

(r\e,jmls) = r" 1 <^(e,r) [y,(fi r ) <g> X s] jm , (18) 

where the wave function <pji{e,r) has usual asymptotic 
at r — > 00. It is either a condition similar to Eq.(T4|) for a 
scattering state or the usual asymptotic with the Whit- 
taker function for a normalizable state (see elsewhere). 
Then using the Wigner-Ekkart theorem, the reduced ma- 
trix element of the C n ^\^ operator may be written as 
follows 

(e K ,jls\\C ntX \\j l s ) = S ssa V n ,x(-) jo+j+lo+r jjo lo 

xC?o°oao W(l s\j;j l) P { $ )x{j0l0) {en, q), (19) 

where quantity V n ,\ is defined by the symmetry prop- 
erties of the spherical harmonics Y\^(yi n ). In particu- 
lar, Vi,x = 1 and V 2 ,x = (-)\ since r x = m " 2 m2 r and 
r 2 = — mi " 1 m2 r, respectively. The function p^ n '(q) reads 

00 

P(ji)\(j„i )( £ K'<l) = J ^i{ e ^ r )3\{qr n )(t)j i {r)dr. (20) 


The quantity C% arising in (|13[) in the case of two-cluster 
system reads 



III. APPLICATION TO DEUTERON ELASTIC 
SCATTERING 

The deuteron is a two-body projectile with only one 
bound state. We analyze the elastic scattering of 2 H as 



a test of the model. Only the s-wave component of the 
ground state wave function is considered. The proton- 
neutron interaction is chosen in a Gaussian form (see 
Table [XJ) with parameters which give the deuteron bind- 
ing energy 2.22 MeV, r.m.s. radius (r 2 ,) 1 ^ 2 = 1.97 fm, 
and triplet scattering length a t — 5.46 fm, close to the 
experimental data. 

As already mentioned, the effect of the deuteron break- 
up on elastic deuteron-nucleus scattering was analyzed 
in many papers [H [H, H [H, [13, [H, S| . We confirm 
here, in particular, that the break-up transition matrix 
elements !F(j i )x(ji) (k, R) fTT]) with even A dominate in 
deuteron induced reactions. Odd partial waves are al- 
most negligible at energy well above the Coulomb barrier 
due to the P 2 ,x coefficient in Eq. (fT9"]) and the similarity of 
the p— A and n— A interactions. The main contributions 
to the polarization potential come from the transitions 
with transfer angular momenta A = 0,2. The deuteron 
break-up energies up to about 15 MeV certainly play a 
leading role in the expansion of the total wave function. 
It allows one to truncate the sum over A in the polar- 
ization potential (|16p at X ma x = 6 and to perform an 
integration over the excitation energy e K up to 40 MeV. 

In spite of the number of applications of the general- 
ized optical model to the deuteron elastic scattering, a 
detailed study of the subject has not been reported. In 
this section we present the results obtained within the ap- 
proach described above to the deuteron elastic scattering 
at energies of 30-50 MeV/u. We consider reactions with 
heavy and light targets separately because of specific fea- 
tures observed in the case of deuteron elastic scattering 
by light ions. 



A. Scattering by heavy targets 

The cluster-target potentials for each reaction were 
obtained from a phcnomcnological optical model anal- 
ysis of appropriate experimental data on the n— A and 
p— A elastic scatterin g at energies around E n<p = Ed/2 
(see discussion in Ref . [16j ) . We use the optical potential 
parametrization given in Ref. [45] in the fitting procedure 
as a starting parameter set. The obtained parameters 
are listed in Table [T] 

The cluster-folding potential U^ 1 ' presents well known 
properties independent of the target ,46] . It is close to 
the sum of the nucleon-target optical potentials. The re- 
sulting radial dependence of the real part has a Woods- 
Saxon shape with a noticeably larger diffuseness parame- 
ter (a ~ 1 fm) and somewhat reduced in magnitude. The 
same is valid for the imaginary part of EA 1 ). The folding 
of the spin-orbit nucleon-target interactions, which are 
usually chosen as a Thomas form (see Table |T|, results in 
a shape close to the Woods-Saxon form with parameters 
similar to the central part. 

The elastic scattering cross section obtained with the 
cluster-folding potential C/W (dashed curves in Fig|TJ) sig- 
nificantly overestimates the experimental data at scat- 
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TABLE I: Effective interaction potentials. 





V a 


H v 


a v 


Wo" 


W D 


Rw 




V so c 


Rso 


aso 


Rc 


n+p d 


62.105 


1.625 




















n + 12 C 


52.25 


2.57 


0.570 




8.05 


2.57 


0.500 


6.2 


2.29 


0.750 




58 Ni 


42.67 


4.53 


0.750 


7.24 


2.586 


4.88 


0.580 


6.2 


4.26 


0.750 




120 Sn 


38.70 


5.77 


0.750 


7.79 


0.375 


6.22 


0.580 


5.5 


5.18 


0.750 




208 Pb 


33.38 


6.93 


0.750 


6.80 




7.48 


0.580 


6.2 


6.52 


0.750 




P + 12 C 


53.29 


2.57 


0.570 




8.05 


2.57 


0.500 


6.2 


2.29 


0.750 


2.75 


5S Ni 


44.92 


4.53 


0.750 


6.10 


2.214 


5.11 


0.534 


6.2 


3.91 


0.750 


4.53 


120 Sn 


48.46 


5.77 


0.750 


6.65 


3.175 


6.51 


0.627 


5.5 


5.18 


0.750 


5.77 


208p b 


48.66 


6.94 


0.750 


8.23 




8.68 


0.580 


6.2 


6.52 


0.750 


6.93 


d+ 12 C 


94.65 


2.24 


0.800 


3.50 


7.40 


3.19 


0.700 








3.29 



"Real part has the Woods-Saxon form V(r) = —Vof(r,R v ,a v ), 
where function f(r,R, a) = (1 + et 1 — -R)/")" 1 . 

''Imaginary part is chosen in the form W(r) = —Wo f(r, R w , a w ) + 
^a-mWD-^fir, Rw,a w ). 

c Spin-orbit interaction has the Thomas form V$o{ r ) = 
2A T (L-s)V SO iJ:/(r, Rso,a so ). 

d The proton-neutron interaction has the Gaussian form V(r) = 
-V e~r 2 /Rl. 
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FIG. 1: Cross sections for deuteron elastic scattering on dif- 
ferent targets. Dashed lines show the cross sections obtained 
with the cluster-folding potential, while solid lines correspond 
to the calculation with the non-local optical potential. The 
experimental data from Ref. 48, 49, 50] are shown by dots. 



tering angles larger than the nuclear-rainbow angle. It 
means that the imaginary part of the folding is too weak. 
This gives us an estimate of the role of the polarization 
potential C/ (2) . 

We calculate the generalized optical potentials + 
for deuteron elastic scattering by 56 Ni, 120 Sn and 
208 Pb targets at energies 80, 85 and 110 MeV, respec- 
tively. The corresponding theoretical cross sections are in 



good agreement with the experimental data (see Fig[JJ). 

The polarization potential U j£(R,R') is illustrated in 
the top panel of Figf2] for the 208 Pb(d,d) reaction and 
J = L = 25. The DPP is symmetric with respect to R 
and R' , therefore we plot the real and imaginary parts of 
DPP in the same figure. The polarization potential is no- 
ticeably non-local, non-monotonic and also L-dependent. 

(2) 

It makes the polarization potential U )£ (R, R ) compli- 
cated for analysis. The L-independent "weighted mean" 
local polarization potential (WLP) was proposed in [ItJ 
as an alternative for the non-local DPP. In the bottom 
panel of Fig[2]the WLP for the same reaction is shown to- 
gether with the module of the corresponding partial wave 
functions. The wave functions ipL=25 calculated with the 
non-local DPP and with the WLP are indistinguishable 
and are shown by the dashed curve in the top part of the 
bottom panel of FigJ^l Note also that the cross sections 
obtained with the WLP and with the initial non-local 
DPP are almost identical. It is illustrated in FigfT] for 
the case of lead target, where the dotted line shows the 
cross section obtained with the WLP. It allows us to an- 
alyze the properties of the WLP instead of the DPP. 

The relative contribution of polarization potential to 
the real part of the optical potential is rather small. In 

(2) 

particular, the value of TteU\ VLP {R) at the minimum 
around R = 8 fm is about -0.7 MeV and amounts to 3% 
of the folding potential i/W, whereas the contribution of 
the imaginary part of the polarization potential is more 
than 30% of the folding potential at this point. For the 
d + 208 Pb collision the total reaction cross section is an 
— 2.9 b at Ed = 110 MeV, where the folding potential 
gives <j r = 2.73 b and the polarization potential gives 

(DPP) ni7K 

<j R =0.17 b. 

Thus, we may conclude that the generalized optical 
potential model provides an adequate description of the 
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FIG. 2: Top panel: Dynamical polarization potential calcu- 
lated for the deuteron elastic scattering on the 208 Pb target 
at Ed =110 MeV and J=L=25. The real (imaginary) compo- 
nent of the DPP is shown by the solid (dash) contours. The 
values of the potential at the extremum are indicated. Bottom 
panel: the module of the partial wave function ipjL=25(k, R) 
corresponding to the d(110 MeV) + 208 Pb reaction is shown in 
top. Solid and dash curves correspond to the wave functions 
calculated with the U {1) cluster-folding and with (U (1) +U (2) ) 
non-local optical potentials. Real and imaginary parts of 
the folding potential (solid curves) and WLP (dash curves) 
is shown for the same reaction in middle and bottom part. 



light two-cluster projectile elastic scattering by heavy nu- 
clei at intermediate energies. The deuteron polarization 
due to the coupling with the break-up channels properly 
describes the missing part of the total reaction cross sec- 
tion (about 10%). 



B. Scattering by light targets 

We study here the d + 12 C collision at 56 MeV, because 
there are experimental data both on the elastic scatter- 
ing [5l| and on the proton-neutron correlations in the 
12 C(d,pn) break-up reaction (52J. We use the parameters 
of the nucleon-carbon potentials taken from Ref . [HJ (see 
TableEJ. 

The general properties of the optical potential are the 




30 60 90 120 150 1 

d cm , deg 



FIG. 3: Elastic scattering cross sections (a) and S-matrices 
(b) for the 12 C(d, d) reaction at 56 MeV. Dash-dotted and 
solid curves in both panels correspond to the results obtained 
with cluster-folding and non-local optical potentials. The long 
dashed line in panel (a) shows the cross section obtained by 
neglecting the polarization potential at L < 6. The short 
dashed line in panel (b) is an S'-matrix generated by phe- 
nomenological optical potential. Dots are the experimental 
data taken from I51I1. 



same as in the case of reactions with heavy targets. The 
elastic scattering cross section calculated with the fold- 
ing potential exceeds the experimental data, see FigO^a). 
However, in contrast with heavy targets an addition of 
the polarization potential to the folding does not lead 
to agreement with experimental data. Comparison of 
the partial S-matrix elements generated by the folding 
(dash-dotted line) and by the non-local optical potential 
(solid line) shows that the polarization potential provides 
a strong additional absorption at low orbital momenta. 

We perform a fit of experimental data on the d + 
12 C elastic scattering within the usual phenomenologi- 
cal optical model [53| using the cluster-folding interac- 
tion as an initial approximation. The fitted OMP pa- 
rameters (see Table HJ provide an angular distribution 
which agrees with the experimental points. The corre- 
sponding S'-matrix elements are shown in Fig|3j6) by 
the short-dashed curve. As it can be seen, there is a 
good agreement of the partial S'-matrix elements obtain- 
ing with the generalized and phenomenological optical 
potentials at L > 6 and significant differences at L < 6. 
Also the phenomenological Sj l elements are rather close 
to those obtained with the folding at low values of L, 
i.e. the break-up probability for central collisions should 
be small. Thus, one may conclude that the model does 
not describe properly the deuteron break-up and, conse- 
quently, elastic scattering at low partial waves. 



Deuteron break-up within the prior-form DWBA 



The function 3~(ji)\(j i ){k, R) in the polarization po- 
tential (fTB)) determines also the prior-form of the DWBA 
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FIG. 4: Left panel: energy spectra of protons in coincidence with neutrons emitted at O n ,lab = 15° for the C(d,pn) elastic 
break-up at 56 MeV in the angular region —60° < 8 P: i a b < 60°. Right panel: angular distribution of p-n correlations in the 
same reaction. Solid and dashed curves are the DWBA calculations of the break-up cross section. Dashed lines show results 
by omitting the contribution of the projectile-target partial waves with L < 6. The calculated triple differential cross sections 
at positive 9 p j a b are renormalized by the factor shown near each curves. Dots are the experimental data from 



break-up amplitude [54 



T 



DW 



(r)|VR, r |^. s .(r)^ +) (R)>, (21) 



where 4> g . s . (r) and 



(r) are the ground and excited 

states of the projectile, and V>p (R) are the distorted 
waves describing the projectile-target relative motion 
in the entrance and exit channels. We calculated the 
deuteron break-up cross section for the 12 C(d,pn) re- 
action at 56 MeV and compared it with experimental 
data [52| , where the angular and energy distributions of 
the protons were measured in coincidence with neutrons 
emitted at 6 n ,iab — 15°. Here we used the same OMP pa- 
rameters as in analysis of d+ 12 C elastic scattering above. 

The results are shown by solid lines in FigfJJ The 
DWBA amplitude (|21[) gives a good agreement with the 
data for negative proton angles and fails at positive ones. 
One may suppose that the differences between the calcu- 
lations and experimental data on deuteron break-up and 
elastic scattering have the same origin. 




FIG. 5: Averaged correlation angle between the protons and 
neutrons emitted in the 12 C(d,pn) break-up reaction at Ed 
= 56 MeV as a function of the impact parameter. The cal- 
culation is performed within the classical three-body model 



2. Deuteron break-up within the classical dynamics model 

To confirm this assumption we performed an analy- 
sis of this reaction within the few-body classical molecu- 
lar dynamics [3]. The models based on the Newtonian 
equations have been successfully applied to the study of 
heavy-ion fragmentation at intermediate energies (see, 
for example, [ll, HE HI, HE HI). Note also that clas- 
sical dynamics approaches turn out to be very effective 
in combination with quantum consideration, allowing to 



explain many aspects of nuclear dynamics using a "tra- 
jectory" language. 

Within the classical model the two-body projectile 
d=(p+n) and target 12 C are treated as classical parti- 
cles moving along the classical trajectories determined by 
the same interactions (Vr jF + Vr) as in the quantum case 
considered above. We tested 10 6 trajectories with ran- 
domly distributed initial parameters (see details in [Oj]). 
Fig[5] shows the averaged angle (A8 pn ) between the pro- 
ton and neutron emitted during the deuteron break-up 
as a function of the impact parameter b — L/p. As it 
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can be seen, protons emitted at negative angles (relative 
to 9 n j a b = 15°, large A0 pn ) are originated mainly in pe- 
ripheral collisions with L > 6, whereas positive angles 
p correspond to central collisions (small (A6 pn ) values). 
Hence, the disagreement with experiment at positive val- 
ues of Op (FigHJ originated from a wrong treatment of the 
contribution of small angular momenta (L < 6) - just as 
in the analysis of elastic scattering. 

Note that a large (A0 pn ) value in peripheral deuteron 
break-up process is caused by the repulsion of protons by 
the Coulomb field of the target, whereas the neutron is 
deflected by attractive nuclear forces. In central collisions 
the effect of the Coulomb forces is much weaker, therefore 
the (A0 pn ) angle turns out to be relatively small. 

The origin of the critical orbital momentum L — 6 (for 
d+ 12 C at 56 MeV) also has a clear explanation in the 
classical model. Using the phenomenological optical po- 
tential for the d(56 MeV) + 12 C reaction (see TableJD, 
and employing an appropriate computational code [591 ] 
we calculated the classical deflection function and sur- 
vival probability 



P S (L) = exp 



W(R)dR 



U(R)] 



(22) 



as a function of the orbital momentum L — bp. P S (L) 
is the probability that the projectile, moving along the 
trajectory with a given impact parameter b, remains in 
the elastic channel. U{R) and W{R) in ((22j) are the 
real and imaginary parts of the optical potential, and 
the integration is performed along the trajectory. The 
calculated survival probability P s are very similar to the 
partial S'-matrix elements, which have the same physi- 
cal meaning (see Figj6]). The deflection function reveals 
the nuclear rainbow angle Onr ~ 70° close to the exper- 
imental value. As it can be seen from Fig[6j the orbital 
momentum Lnr = 6 corresponds to the nuclear rain- 
bow scattering. This means that the trajectories with 
L < L pass deeply in the interaction region, while the 
trajectories with L > Lnr are more peripheral. Thus 
the model used here for a calculation of the generalized 
optical potential does not treat properly central collisions 
with a strong overlapping of the colliding nuclei. 



3. Projectile-target non-physical bound states 

The phenomenological OM analysis shows that the ab- 
sorption at low orbital momenta in the elastic scattering 
is well described by the imaginary part of the folding 
potential only and does not require any addition (com- 
pare dash-dotted and dashed curves in Fig[3j6)). Note 
that a simple cut of the polarization potential at L < 6 
in the elastic scattering analysis and omitting the con- 
tribution of these partial waves to the DWBA break-up 
cross-section lead to a significant improvement of the re- 
sults in both cases. It is shown by the dashed curves in 




FIG. 6: Deuteron survival probability and deflection function 
calculated within classical model [591 ] for the d(56 MeV) + 12 C 
reaction. Open squares show the Sjl matrix obtained within 
the phenomenological optical model with the same potential. 



FigGUa) and Fig0J Thus, the experimental data indi- 
cates that the deuteron in the inner region of the target 
nucleus turns out to be stable relatively to the break-up 
channels, that is confirmed also in previous studies [25| . 

Thus, we conclude that the calculated DPP overes- 
timates the absorption at small partial waves, i.e. at 
small projectile-target relative distance. The radial de- 
pendence of the polarization potential in this region is 
defined mainly by the partial Green function g^£ (R, R') 
(fTS"]) . The properties of the Green function at low L val- 
ues are significantly affected by the properties of the fold- 
ing potential U^' (R) which is used to calculate the Green 
function. The interaction (R) is calculated as a sum 
of the folded complex cluster-target optical potentials. 
Parameters of these potentials are usually fitted in order 
to reproduce experimental cross sections. In this pro- 
cedure the scattering phase-shifts are retrieved but not 
the wave functions. Thus, the cluster-target potentials 
may provide an incorrect behavior of the partial wave 
functions at small distances because of ir ambiguity of 
the phases. In particular, cluster-target potentials and, 
consequently, cluster-folding one turn out to be deep and 
contain many forbidden bound states. This may results 
in incorrect radial dependence of the Green function at 
small distance, since the corresponding partial wave func- 
tions penetrate deeply into the interaction region. This 
leads to the rise of the DPP at small R values. 

The observed stability of deuteron moving in nuclear 
matter with respect to the break-up means in fact that 
the deuteron does not penetrate deeply into a target due 
to the Pauli blocking. The effects of antisymmetrization 
in deuteron elastic scattering have been studied before 
[6fj| . In order to take it into account consistently within 
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FIG. 7: (a) Imaginary part of the partial Green functions (L—2) and real part of corresponding interaction potentials for 
the 12 C(d,d) reaction. Functions gL and g*j? are calculated using effective potential Vl = + ti 2 L(L + l)/2mR 2 and 
its supersymmetric partner V^ S ' after removal of the state E2 = —(18.5 + z9.3) MeV, respectively, (b) WLPs for the same 
reaction resulted from calculation with non- modified non-local optical potential (dotted line), with non-local potential neglected 
polarization part at L < 6 (dash line), and with non-local optical potential after removal forbidden states (solid line). 



our approach one needs to remove the forbidden states 
from the nucleon-target potentials, that makes them non- 
local and results in the complication of their further 
treatment. Therefore we apply a simplified method mod- 
ifying the d— 12 C cluster-folding potential U^(R) which 
also presents a number of non-physical bound states. Let 
us then remove these states. 

States found in the d— 12 C folding potential (without 
spin-orbit interaction) are listed in Table HT1 E' nL are the 
eigenvalues corresponding to the states in the potential 
without imaginary part, while E n L are the eigenvalues 
in the complex potential (since the cluster-folding inter- 
action is complex). The imaginary part of the potential 
leads to the appearance of a negative imaginary addition 
to the eigenvalues. He E n L of the bound states as well 
as of the narrow resonances are modified a little, while 
the broad resonances are shifted significantly. 

For an hermitian Hamiltonian the S'-matrix poles cor- 
responding to resonant states are symmetric with respect 
to the imaginary p-axis. In the case of a complex po- 
tential this symmetry is broken. Generally the resonant 
poles move in clockwise direction in the complex p-plane 
(see Ref.[6l| for details). The states with IZep > 
(right half-plane) get negative addition to the Imp, while 
some of the poles in the left complex half-plane cross the 
real axis and become normalizable states (IZep < and 
Imp > 0, i.e. iPl(R — > 00) ~ e" Im|,fl ). The energies 
corresponding to these states in the d— 12 C folding po- 
tential are also listed in Table HT1 

We apply the technique explained in [62|, [63| , which al- 
lows to eliminate normalizable states from the spectrum 
of the complex potential using supersymmetric trans- 
forms in each partial wave. The resulting potential be- 
comes L-dependent, and contains a strong repulsive core 



TABLE II: Bound, resonant and normalizable states in the 
d + 12 C folding potential. 



L 


nodes 


K L (MeV) 


E nL (MeV) 








-51.47 + iO 


-51.29 - «8.25 


1 





-34.60 + iO 


-34.52 - i9.ll 





1 


-19.72 + iO 


-19.57 - i8.34 


2 





-18.57 + iO 


-18.48 - i9.25 


1 


1 


-6.93 + iO 


-6.37 - i6.82 


3 





-4.03 + iO 


-3.62 - i8.39 





2 


0.14 - tO.O 


0.49 - i0.65 


1 


2 


1.24 - i2.10 


0.82 - il.25 


2 


1 


1.89 - iO.18 


1.76 - il.72 


3 


1 


5.95 - i4.51 


3.51 - i3.49 


4 





7.36 - i0.59 


7.11 - i3.57 


5 





16.45 - i4.88 


12.12 - i7.45 


0" 






1.47 - i3.06 


2 a 






3.85 - i3.57 


4" 






9.57 - i5.78 


"The 


normalizable 


states in the d + 12 C foldinj 


; potential with 



positive real energy. 



at small distances. The supersymmetric transforms do 
not modify the scattering phase-shifts. Nevertheless the 
partial wave function and, consequently, the Green func- 
tion turns out to be pushed out from the interaction re- 
gion (see FigJT^a)). The obtained Green function was 
used in the DPP calculation. The corresponding WLP 
are shown in FigJTf 6) in comparison with the initial WLP 
and with the WLP resulted from a dropping of the non- 
local polarization potential at L < 6. We may conclude 
that the elimination of the forbidden states modifies the 
DPP in a correct way. 

Modules of the 5-matrix elements for the same reaction 
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FIG. 8: The S-matrices for 1 C(d,d) reaction at energy Ed = 56 MeV are compared in panel (a). Full and open squares 
show the S-matrix generated by non-modified non-local and empirical optical potential, respectively. Full and open circles are 
the S-matrix elements calculated with non-local optical potential after removal forbidden states and with the corresponding 
WLP, respectively, (b) Elastic scattering cross section for the same reaction. The solid curve shows the calculations with the 
WLP after removal of the forbidden states, while the dashed curve shows the calculation with non-modified non-local optical 
potential. Dots are the experimental data [511 ]. 



are shown in Fig[5Ja). One may see that the supersym- 
metric transform leads to a decrease of the absorption 
in elastic channels at low partial waves. Thus the elim- 
ination of the non-physical states in the d— 12 C folding 
potential allows one to describe effectively the suppres- 
sion of the deuteron break-up at low values of angular 
momenta. 

The d— 12 C folding interaction does not contain nor- 
malizable states with L — 5 as it may be expected 
from the behavior of the phcnomcnological S'-matrix (see 
FigJS^a)). This indicates that the folding potential is not 
the best substitution for the VR.r hi the Green function 
calculation. Damped Sl=5 matrix element (solid circles 
in FiglUJa)) results in the oscillating behavior of the cross 
section at large scattering angles. This problem is over- 
came somehow if we use the WLP (solid line in Fig[7j6)), 
which smoothes the S'-matrix (open circles in FigJSJJa)) 
by averaging the non-local polarization potential over all 
orbital momenta. The angular distribution in the d+ 12 C 
elastic scattering is shown in FigJ5]^6) together with ex- 
perimental data and with the cross section obtained with- 
out the supersymmetric transforms. 



IV. CONCLUSIONS 

By extending the model proposed in earlier papers, we 
derive the generalized optical potential for elastic scat- 
tering of a few-cluster projectile, taking into account ex- 
plicitly the coupling with the break-up channels. We do 
not use most of the simplifications which were employed 
in previous papers. In particular, applying the model 
to deuteron elastic scattering we take into account the 



spin of projectile, consider the coupling to the projectile 
continuum with cluster-cluster relative orbital momenta 
I < 6, and apply the suitable approximation of the Green 
function instead of the free-particle one used before. It al- 
lows to improve an agrement with the data and previous 
results obtained within different approaches, that sup- 
ports the efficiency of the model. The model was applied 
to study of the deuteron elastic scattering at energies of 
few tens of MeV per nucleon and good agreement with 
experimental data was obtained. 

It was also shown that the behavior of the polarization 
potential at low orbital momenta is noticeably affected 
by the non-physical bound states in the projectile-target 
system. In the case of light heavy-ion scattering it leads 
to the overestimation of the absorption in the GOP in 
low partial waves. The elimination of these forbidden 
states allows one to obtain an appropriate polarization 
potential. 

Note that the non-physical bound states do not reveal 
itself in the deuteron scattering by heavy targets in spite 
of their existence for low partial waves. The reason is 
a much stronger absorption part of the folding potential 
in the case of heavy targets as compared with light nu- 
clei. Addition of the polarization potential to the folding 
one gives a negligible effect on the elastic scattering cross 
section in low partial waves because the contribution of 
these partial waves is suppressed by the absorptive part 
of the folding interaction. Note, however, that the forbid- 
den states may play some role in other reaction channels, 
for example, in break-up. This subject is an interesting 
problem for future studies. 

Application of the proposed model to the reactions 
with a three-cluster weakly bound nuclei (such as 6 He 
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= a + n + n) will be done in future works. 

Acknowledgments 

This work was supported by the INTAS (Grant No. 
04-83-2649) and by a grant of FNRS (Belgium). One 



of the authors (A.D.) would like to thank the PNTPM 
department of the Universite Libre de Bruxelles for hos- 
pitality. Dr. J.-M. Sparenberg is kindly appreciated for 
useful discussions and providing the SUSY transforma- 
tion code. 



B. Jonson, Phys. Rep. 389 (2004) 1 [31 
N. Keeley, R. Raabe, N. Alamanos, J.L. Sida, Prog. Part. 
Nucl. Phys. 59 (2007) 579 

Y. Suzuki, R.G. Lovas, K. Yabana, K. Varga, Structure [32 
and reactions of light exotic nuclei, Taylor and Francis, 
2003 [33 
M.V. Zhukov, B.V. Danilin, D.V. Fedorov, J.M. Bang, [34 
I.J. Thompson and J.S. Vaagen, Phys. Rep. 231 (1993) 
151 [35 
B.V. Danilin, I.J. Thompson, J.S. Vaagen, M.V. Zhukov, 
Nucl. Phys. A 632 (1998) 383 [36 
B. V. Danilin, J. S. Vaagen, T. Rogde, S. N. Ershov, I. J. 
Thompson, and M. V. Zhukov, Phys. Rev. C 73 (2006) [37 
054002 

P. Descouvemont, E. Tursunov, D. Baye, Nucl. Phys. [38 
A765 (2006) 370 [39 
S.N. Ershov, B.V. Danilin, and J.S. Vaagen, Phys. Rev. 
C 64 (2001) 064609 [40 
Yu. Ts. Oganessian, V. I. Zagrebaev, J. S. Vaagen, Phys. [41 
Rev. C 60 (1999) 044605 [42 
R. Chatterjee, P. Banerjee, R. Shyama, Nucl. Phys. 
A692 (2001) 476 [43 
K. Rusek, K. W. Kemper, R. Wolski, Phys. Rev. C 64 [44 

(2001) 044602 

R.S. Mackintosh, N. Keeley, Phys. Rev. C 70 (2004) [45 
024604 

B. Abu-Ibrahim and Y. Suzuki, Phys. Rev. C 70 (2004) [46 
011603(R) [47 
V.I. Zagrebaev, D.N. Semkin, Bull. Rus. Acad. Sci. Phys. 
59, No.5 (1995) 145; ibid. 59, No.ll (1995) 140 [48 
R.G. Johnson, P.J.R. Sopper, Phys. Rev. C 1 (1970) 976 [49 

G. Rawitscher, Phys. Rev. C 9 (1974) 2210 [50 
Y. Sakuragi, Phys. Rev. C 35 (1987) 2161; Y.Sakuragi [51 
et al., Prog. Theor. Phys. 70 (1983) 1047 [52 
N. Austern, Y. Iseri, M. Kamimura, M. Kawai, G. Raw- [53 
itscher and M. Yahiro, Phys. Rep. 154 (1987) 125 
K. Hagino, N. Rowley and AT. Kruppa, Comput. Phys. [54 
Commun. 123 (1999) 143 [55 
V.I. Zagrebaev and V.V. Samarin, Yad. Phys. 67 (2004) 

1488; Phys. At. Nucl. 67 (2004) 1462 [56 
A.Diaz-Torres and I.J. Thompson, Phys. Rev. C 65 [57 

(2002) 024606 

T. Matsumoto et al., Phys. Rev. C 70 (2004) 061601(R) [58 
T. Matsumoto et al., Phys. Rev. C 73 (2006) 051602(R) 

H. Feshbach, Ann. Phys. 5 (1958) 357; ibid. 281 (2000) [59 
519 

J. Testoni, L.C. Gomes, Nucl. Phys. 89 (1966) 288 [60 

G. Baumgartner, Z. Phys. 204 (1967) 17 

V.I. Zagrebaev, Izv. AN SSSR, Ser. Fiz. 42 (1978) 131 

S. Mukherjee, Nucl. Phys. A118 (1968) 423 

W.G. Love et al., Nucl. Phys. A291 (1977) 183 [61] 

M. Hussein et al., Phys. Rep. 113 (1984) 133 



A. Bouyssy, N. Vinh Mau, D.M. Brink, Phys. Lett. B102 

(1981) 109; N. Vinh Mau et al., Nucl. Phys. A457 (1986) 
413 

H. Dermawan and F. Osterfeld, V. A. Madsen, Phys. 
Rev. C 25 (1982) 180 

J.C. Pacheco et al., Nucl. Phys. A588 (1995) 537 
P.W. Coulter, G.R. Satchler, Nucl. Phys. A293 (1977) 
269 

S.B. Khadkikar, L. Rikus, A. Faessler, R. Sartor, Nucl. 
Phys. A369 (1981) 495 

M.S. Hussein, MP. Pato, L.F. Canto, R. Donangelo, 
Phys. Rev. C 47 (1993) 2398 

G. R. Satchler, Direct Nuclear Reactions, Capther 4.3, 
Oxford University Press, 1983 

S. Watanabe, Nucl. Phys. A8 (1958) 484 

H. Nishioka, J. A. Tostevin, R.C. Johonson, K.-I. Kubo, 
Nucl. Phys. A415 (1984) 230 

F. Petrovich et al., Nucl. Phys. A425 (1984) 609 

H. Amakawa, K.-I. Kubo, Nucl. Phys. A266 (1976) 521 

R.S. Mackintosh and A.M. Kobos, Phys. Lett. B116 

(1982) 95 

H. Amakawa and K. Yazaki, Phys. Lett. B87 (1979) 159 
A. Faessler, L. Rikus, R. Sartor, Comp. Phys. Comm. 28 

(1983) 275 

CM. Perey, F.G. Perey, At. Data Nucl. Data Tab. 13 
(1974) 293 

F.G. Perey, G.R. Satchler, Nucl. Phys. A97 (1967) 515 

I. J. Thompson, M.A. Nagarajan, J.S. Lilley, M.J. Smith- 
son, Nucl. Phys. A505 (1989) 84 

E. J. Stephenson et al., Phys. Rev. C 28 (1983) 134 
J. Bojowald et al., Phys. Rev. C 38 (1988) 1153 
A.C. Betker et al., Phys. Rev C 48 (1993) 2085 

N. Matsuoka, et al., Nucl. Phys. A455 (1986) 413 
N. Matsuoka, et al., Nucl. Phys. A391 (1982) 357 
V.I. Zagrebaev, A.S. Denikin, A. P. Alekseev, Optical 
model code of NRV, |http: / /nrv.jinr.ru/nrv| 

F. Rybicki and N. Austern, Phys. Rev. C 6 (1972) 1525 

G. Royer, Y. Raffray, A. Oubahadou, and B. Remaud, 
Nucl. Phys. A466 (1987) 139 

R. Wada et al., Nucl. Phys. A539 (1992) 316 

K. Mohring, T. Srokowski, D.H.E. Gross, Nucl. Phys. 

A533 (1991) 333 

A.S. Denikin, V.I. Zagrebaev, Yad. Fiz. 65 No.8 (2002) 
1494; Phys. At. Nucl. 65 No.8 (2002) 1459 
V.I. Zagrebaev, A.S. Denikin, A. P. Alekseev, Classical 



model code of NRV, http://nrv.jinr.ru/nrv 
W.S. Pong, N. Austern, Ann. Phys. (NY) 93 (1975) 369; 
A. A. Ioannidis, R.C. Johnson, Phys. Rev. C 17 (1978) 
1331; J.A. Tostevin, M.H. Lopes, R.C. Johnson, Nucl. 
Phys. A465 (1987) 83 

J.-M. Sparenberg and D. Baye, Phys. Rev. C 54 (1996) 
1309 



13 



[62] D. Baye, Phys. Rev. Lett. 58 (1987) 2738; J. Phys. A20 (1996) 435 

(1987) 5529 

[63] D. Baye, G. Levai, J.-M. Sparenberg, Nucl. Phys. A599 



